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Abstract 

Semiclassical approximations for tunneling processes usually involve complex trajectories or com- 
plex times. In this paper we use a previously derived approximation involving only real trajectories 
propagating in real time to describe the scattering of a Gaussian wavepacket by a finite square 
potential barrier. We show that the approximation describes both tunneling and interferences very 
accurately in the limit of small Planck's constant. We use these results to estimate the tunneling 
time of the wavepacket and find that, for high energies, the barrier slows down the wavepacket but 
that it speeds it up at energies comparable to the barrier height. 

PACS numbers: 03.65.Sq,31.15.Gy 
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I. INTRODUCTION 



The success of semiclassical approximations in molecular and atomic physics or theoret- 
ical chemistry is largely due to its capacity to reconcile the advantages of classical physics 
and quantum mechanics. It manages to retain important features which escape the classical 
methods, such as interference and tunneling, while providing an intuitive approach to quan- 
tum mechanical problems whose exact solution could be very difficult to find. Moreover, 
the study of semiclassical limit of quantum mechanics has a theoretical interest of its own, 
shedding light into the fuzzy boundary between the classical and quantum perspectives. 

In this paper we will apply the semiclassical formalism to study the scattering of a one 
dimensional wavepacket by a finite potential barrier. In the case of plane waves, the tunneling 
and refiection coefficients can be easily calculated in the semiclassical limit, giving the well 
known WKB expressions For wavepackets, however, the problem is more complicated 



and few works have addressed the question from a dynamical point of view 



3|, IJ] . The time 



evolution of a general wavefunction with initial condition tp{x, 0) = ipo{x) can be written as 

ij{x,T) =< x\K{T)\tlJo >= J < x\K{T)\xi > dxi < a;#o >, (1) 

where K{T) = e~*^^/^ is the evolution operator and H is the (time independent) hamil- 
tonian. The extra integration on the second equality reveals the Feynman propagator 
< x\K{T)\xi >, whose semiclassical limit is known as the Van-Vleck formula Q (see 
next section). When the Van-Vleck propagator is inserted in Eq.([T]) we obtain a general 
semiclassical formula which involves the integration over the 'initial points' xf. 



1psc{x,T) = j < x\K{T)\Xi >VanVlcck dXi < 



'0 > . (2) 



If this integral is performed numerically one obtains very good results, specially as h goes 
to zero. However, doing the integral is more complicated than it might look, because for 
each Xi one has to compute a full classical trajectory that starts at Xi and ends at x after a 
time T, which may not be simple task. Alternative methods involving integrals over initial 



conditions (instead of initial and fina 



and shown to be very accurate aaa. 

All these approaches sum an infinite number of 
contributions and hide the important information of what classical trajectories really matter 
for the process. 



coordinates) in phase space have also been developed 
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In a previous paper |9| several further approximations for this integral were derived 
and applied to a number of problems such as the free particle, the hard wall, the quartic 
oscillator and the scattering by an attractive potential. The most accurate (and also the most 
complicated) of these approximations involves complex trajectories and was first obtained 



by Heller and collaborators [10|, lUl- The least accurate (and the simplest to implement) 
is known as the Frozen Gaussian Approximation (FGA), and was also obtained by Heller 
[l2|. It involves a single classical trajectory starting from the center of the wavepacket. 



These 



However, other approximations involving real trajectories can be obtained 
are usually not as accurate as the complex trajectory formula, but are much better than 
the FGA and can be very good in several situations. Moreover, it singles out real classical 
trajectories from the infinite set in Eq.([2]) that can be directly interpreted as contributing 
to the propagation. 

In this paper we apply these real trajectory approximations to study the tunnel effect. 
Since this is a purely quantum phenomena, it is a very interesting case to test the semi- 
classical approximation and to understand what are the real trajectories that contribute 
when the wavepacket is moving 'inside' the barrier. More specifically, we will consider the 
propagation of a Gaussian wavepacket through a finite square barrier. We shall see that the 
semiclassical results are very accurate, although some important features of the wavepacket 
propagation cannot be completely described. 

This paper is organized as follows: in the next section we review the semiclassical results 
tl 

derived in [9|, which are the starting point of this work. Next we describe the evolution of 
a Gaussian through a square potential barrier in its three separate regions: before, inside 
and after the barrier. Finally in section IV we discuss the calculation of tunneling times, as 

n 

proposed in [^]. We find that the barrier slows down the wavepacket at high energies, but 
that it speeds it up at energies comparable to the barrier height. Finally, in section V we 
present our conclusions. 

II. APPROXIMATION WITH COMPLEX AND REAL TRAJECTORIES 

One important class of initial wavefunctions is that of coherent states, which are minimum 
uncertainty Gaussian wavepackets. In this paper we shall consider the initial wavepacket 
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IV'o > as the coherent state of a harmonic oscillator of mass m and frequency uj defined by 



1^) =e-^l^l'e^^'|0), (3) 

where |0) is the harmonic oscillator ground state, is the creation operator and z is the 
complex eigenvalue of the annihilation operator a with respect to the eigenfunction 1^;). 
Using the position and momentum operators, q and p respectively, we can write 

1 



■ P 

i- 

c 



V2 \b \ 



(4) 



where q and p are real numbers. The parameters b = (h/muj)^ and c = (hmuj)^ are the 
position and momentum scales respectively, and their product is h. 

In order to write the Van Vleck formula of the Feynman propagator, we need to introduce 
the tangent matrix. Let S = S{xf,T;Xi,0) be the action of a classical trajectory in the 
phase space (X, P), with Xi = X{0) and xj = X{T). A small initial displacement {6xi,Spi) 
modifies the whole trajectory and leads to another displacement {6xf,6pf) at time T. In 
the linearized approximation, the tangent matrix 071 connects these two vectors of the phase 
space 

^ 1 \ /5x,\ I 
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(5) 



where Sa = d'^S/dxf, Sif = d'^S/dxidxf = Sf^ and Sff = d'^S/dx'j. In terms of the 
coefficients of the tangent matrix, the Van Vleck propagator is jsl 

1 



{X f\K{T)\Xi)vanVleck 



exp 



^S{xf,T;Xi,0) - 



(6) 



b^y2nmqp 

For short times rrigp is positive and the square root is well defined. For longer times rugp 
may become negative by going through zero. At these 'focal points' the Van Vleck formula 
diverges. However, sufficiently away from these points the approximation becomes good 
again, as long as one replaces rrigp by its modulus and subtracts a phase tc/2 for every focus 
encountered along the trajectory. We shall not write these so-called Morse phases explicitly. 

Assuming some converging conditions, the stationary phase approximation allows us to 
perform the integral over Xi in Eq. ([2]) (for more details, see j^). We obtain 



ip{z,Xf,T), 



6-1/27,-1/4 

^/niqq + iUlqp 



exp 



'-S{xf, T- xo, 0) + '-p{x, - q/2) - 



(7) 



where xq is the value of the initial coordinate Xi when the phase of the propagator is sta- 
tionary. It is given by the relation 

The end point of the trajectory is still given by X{T) = Xf. In spite of q and p being real, xo 
and Po are usually complex. This implies that the classical trajectories with initial position 
Xo and momentum pq are complex as well, even with xj G M. Eq. ([7]) was first obtained by 



Heller 
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111 ] and it is not an initial value representation (IVR). There are a priori several 
complex trajectories satisfying the boundary conditions. Thanks to the stationary phase 
approximation, we were able to replace an integral over a continuum of real trajectories ([2]) 
by a finite number of complex ones The problem is now solvable, but still quite difficult 
to compute. However, it turns out that, in many situations, these complex trajectories can 



be replaced by real ones, which are much easier to calculate 

Therefore, we look for real trajectories that are as close as possible to the complex ones. 
Let {X(t), P(t)) G C X C be the coordinates of a complex trajectory, and {u(t),v{t)) a new 
set of variables defined by 

1 fx P\ I fx P\ 

According to Eq. ([H]), the boundary conditions become 

-(°' - ;^ (y + 't) - (f -'D - ^(^'-^" 

The initial condition is then the complex coordinate z and the final condition is the real 
position Xf. The real and imaginary parts of z are related to the central position q and 
the central momentum p of the initial wavepacket respectively. This gives us three real 
parameters that we may use as boundary conditions to determine the real trajectory. But 
a particle whose initial conditions are q and p will not a priori reach xj after a time T. 
Although it is possible to satisfy such final condition, it will not usually happen because 
X(T) is imposed by q and p. Likewise, fixing the initial and final positions q and Xf will 
not generally lead to -P(O) = p. Therefore we need to choose only two boundary conditions 
among the three parameters, and use the hamiltonian of the system to calculate analytically 
or numerically the third one. This means that the relation ([8]) will not be generally fulfilled 
and the hope is that it will be fulfilled approximately. For a discussion about the validity of 



this approximation, see the beginning of the third section in 9|. If we fix (g, p), we obtain the 



Frozen Gaussian Approximation of Heller 



12|. This is an initial value representation that 



involves a single trajectory and is unable to describe interferences or tunneling, which are 
the aim of this paper. However, we can fix X{0) = q, X{T) = Xf and calculate P(0) = Pi. 
When the complex quantities in Eq.([7]) are expanded about this real trajectory we obtain 







5-1/2^-1/4 



exp 



im. 



qp 



TS{xf,T;q,0) + ^pq-^ 
n 2h 2 mqq + iniqp 



c 



Eq. (ITT!) is the semiclassical formula we are going to use in this paper. We shall show 
that, although still very simple, it can describe tunneling and interferences quite well. 



III. THE 1-D SQUARE BARRIER 



Consider the specific case of a particle of unit mass scattered by the 1-D square barrier 
defined by (see figH]) 



IVn if X E [—a, a] where a G 
(12) 
otherwise 

The initial state of the particle is a coherent state iplz, x, 0) = {x\z) with average position 
q < —a and average momentum p > 0, i.e., the wavepacket is at the left of the barrier and 
moves to the right. In all our numerical calculations we have fixed Vq = 0.5 and defined the 
critical momentum p = \/2Vq = 1. 

The application of the semiclassical formula EqUT] requires the calculation of classical 
trajectories from qtoxj in the time T. For the case of a potential barrier, the number of 
such trajectories depends on the final position xj. This dependence, in turn, causes certain 
discontinuities in the semiclassical wavefunction. 

Since the initial wavepacket starts from q < —a, it is clear that for > a (at the right 
side of the barrier) there is only one trajectory satisfying x(0) = q and x{T) = xj. This 
'direct trajectory' has pi > a/2Vo and x(t) increases monotonically from q to Xf. 

For Xf < —a, on the other hand, in addition to the direct trajectory there might also be 
a refiected trajectory, that passes through xj, bounces off the barrier and returns to in 
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the time T. The initial momentum of such a reflected trajectory must be greater than that 
of the direct one, since it has to travel a larger distance. However, if this distance is too big, 
i.e., if << —a, the initial momentum needed to traverse the distance in the fixed time 
T becomes larger than \/2Vq and the refiected trajectory suddenly ceases to exist (see next 
subsection for explicit details for the case of the square barrier and figure 2 for examples). 

This qualitative discussion shows that refiected trajectories exist only if is sufficiently 
close to the barrier. The points where these trajectories suddenly disappear represent discon- 
tinuities of the semiclassical calculation. Fortunatelly, this drawback of the approximation 
becomes less critical as h goes to zero, since the contribution of the refiected trajectory at 
those points become exponentially small as compared to the direct one (see for instance 
figure 2(g)). 

In the remaining of this paper we are going to obtain explicit expressions for ip{z, Xf, T)sc 
before, inside and after the barrier. For fixed q we will calculate the classical trajectories for 
each Xf, extracting the initial momentum pi, the action S = S{xf, T; q, 0) and its derivatives 
(in order to obtain mqg and niqp). 



A. Before the barrier: Xf < —a 



The specificity of this region is that there may exist two different paths connecting q to 
X during the time T: a direct trajectory and a reflected one (flg. [1]) whose initial momenta, 
action and tangent matrix elements are given by 

{x -qf . T 



Pi 



Pid — Tj, 1 

X + q + 2a 



' ■99 

\2 



m„„d = 1; m. 



qpc 



A 



{x + q + 2ay 



- 1 , Tflqp r 



T 
A' 



(13) 
(14) 



T 2T 
where X = b/c. The contribution of each trajectory to the wavefunction at Xf, ipd and ipr, is 



^rf = — 7^=^ exp 



i jxf - qf 
h 2T 



1 iT f pT — Xf + q 



2h^^ 2 X + iT 



cT 



iPr = / .r exp 



^ i(xf + q + 2af i 1 iT fpT + Xf + q + 2a 
'^^h 2T + 2^^^-2AT^i ^ 



(15) 

Notice that we have added an extra phase Omipr. Without this extra phase (that includes 
the minus sign coming from the tangent matrix elements in Eq. (ll4l) ). the wavepacket would 



not be continuous as it goes through the barrier. For a hard wall, for instance, we impose 
6' = TT to guarantee that the wavefunction is zero at the wall. For smooth barriers this 
phase would come out of the approximation automatically, but for discontinuous potentials 
we need to add it by hand. To calculate 9 we rewrite the previous expressions in complex 
polar representation, V'd = D{xf)e'^''^''f\ ipr = i?(x/)e*'^'-(^/)+*^ and let W'(a;/)e*^'-(^/)+^« be 
the wavefunction inside the barrier, where ^ is the corresponding phase correction. The 
continuity of the wavefunction at x/ = —a imposes 



(16) 



Eqs. f fTSj) show that R{—a) = D{—a) and ipdi—a) = </?,,(— a). Denoting ip = ipyj{—a) — 
<^d{-a), Eq. (HI]) becomes 1 + e*^ = v4e*(^+^) where A = [W{-a)/ D{-a)]. This complex 
relation represents in fact two real equations for the unknown variables 6 and ^. The solutions 
consistent with the boundary conditions are cos(^) = — 1 and cos(0 + ^) = A/2. In 
the limit where p goes to zero (or the potential height Vq goes to infinity) we obtain 6 = n 
as expected. Finally, the full wavefunction before the barrier is ilj{z,Xf,T)sc = ipd + i^r and 
the probability density can be written as 



[z,Xf,T), 



1 



1 



1 + 



rp2 



exp 



A2 + r2 



Xf — q — pT 



+ exp 



+ 2 cos 



A^ 



A2 + T2 

2{xf + a) 
h{X^ + T2) 



Xf + q + pT + 2a 



{X^p - (g + a)T) - 9 



exp 



A^ 



{pT + q + af + {Xf + a 



i2n 



A2 + T2 



62 



(17) 



This is the same result as obtained in [9] for a completely repulsive barrier (Vq — > oo), 
except for the phase, because of the different boundary condition aXxf = —a {\ip{—a)sc\'^ = 
for the hard wall). However, as discussed in the beginning of this section, an additional 
difficulty appears when the wall is finite: the refiected trajectory does not always exist. 
From the classical point of view, there is no refiected part if the energy E = pfr/2 > Vq. The 
maximum initial momentum allowed is then ^/2Vo and a particle with such momentum takes 



the time Tr 



a + q 
7^ 



to reach the barrier. Furthermore, for T > Tc the refiected trajectory 



only exists if pi 



Xf + q + 2a 



^ i.e. if \xf\ 



— Xf ^ Xr 



q + 2a + y/WoT. 



Therefore, if T ^ Tc and ^ Xc, the probabihty density is given by Eq. (fT7|) . otherwise we 
only have the contribution of the direct tpd and 

/ Xf — q — pT^ 



',Xf,T)g(.\ — — exp 



A2 + T2 



(18) 



As a final remark we note that the calculation of 6 might involve a technical difficulty 
depending on the values of h, p and T. For some values of these parameters the contribution 
of the direct and reflected trajectories might become very small at x=-a (see for instance 
flg.2(f), which shows the reflected wavepacket in a case of large transmission). In these 
cases the probability density becomes very small at x=-a and the value of the phase 6 is 
irrelevant. In some of these situations, where the value of 6 does not affect the results, we 
actually found that cos{theta) = — 1 > 1, which cannot be solved for real 6. For the 
sake of numerical calculations we have set ^ = in these cases. 

The semiclassical wavepacket is now completely described for Xf < —a. The probability 
density lipsd"^ is a function of q,p,Xf,T and depends on several parameters, a,b, h and Vq. 
In our numerical calculations we flxed a = 50. This makes the barrier large enough so 
that we study in detail what is happening inside (see subsection IIII Bl) . The height of the 
barrier intervenes only in Tc and Xc to establish the limits of the reflected trajectory. Its 
numerical value is not important, but its comparison with p is fundamental: since we have 



fixed Vq = 0.5 this gives Pir ^ p = V2Vq = 1. Finally, to simplify matters we fixed b = c, 
i.e. the same scale for position and momentum. This imposes X = b/c = 1. Quantum 
phenomena such as interference and tunneling should be more important for high values of 
h. Since h = be = b^, b becomes in fact the only free parameter of the approximation. We 
have also fixed q = —60, which guarantees that the initial wavepacket is completely outside 
the barrier for all values of b used. 

Fig. [2] shows snapshots of the wavepacket as a function of Xf at time T = 50. Consider 
first the panels (a)-(c) with h = 1. The agreement between the exact and the semiclassical 



curves is qualitatively good for p ^ V2Vo = 1. The interference peaks occurs at about the 
same positions, but the height of the peaks are not exactly the same. Also the intervals 
between peaks are a little bigger for the semiclassical curve than for the exact one. On the 
other hand, when p is increased, the comparison gets worst and the approximation is not 
really accurate for p = 2. However, we see that the value of iV'scP at p = 2 is only a tenth 
of its value at p = 0.5: the most important part of the wavepacket is in fact inside and after 



the barrier. It is then really important to consider xj > —a for high p and we need to wait 
until subsections IIIIBI and IIII CI to look at the whole picture. 

When h = 0.25, Fig. [2](d)-(f) and (h), the approximation improves substantially, espe- 
cially close to the barrier; this shows that the extra phase 6 works well. When p is increased, 
the contribution of the direct trajectory becomes irrelevant and the interference oscillations 
are lost in the semiclassical calculation, although it still shows good qualitative agreement 
in the average. The cut-off of the semiclassical curve at x/ = —Xc is also clearly visible, 
whereas the exact one is decreasing continuously. On the one hand this means that the 
approximation is not perfect but, on the other hand, the semiclassical approximation ex- 
plains that the fast rundown of the exact quantum wavepacket comes from the progressive 
disappearance of the reflected classical trajectory due to the finite size of the barrier. Fi- 
nally, for h = 0.1, Fig. El^g), the approximation becomes nearly perfect. As expected, the 
semiclassical approximation works better and better when h is decreasing, i.e. when the 
quantum rules give way to the classical ones. 

To end this subsection, we mention that the quality of the approximation is independent 
of the time T, except for times slightly smaller than Tc. In this time interval only the direct 
trajectory contributes but the exact wavepacket already shows interferences that can not be 
described by {ipscl"^ (fig- El T = 10). We now enter the heart of the matter, and consider 
what's happening inside and after the barrier. 

B. Inside the barrier: —a ^ xj ^ a 

From the classical point of view there is only the direct trajectory in this region (see Fig. 
[1]), since a reflection on the other side of the barrier (at x = a) can not be considered without 
quantum mechanics. Calling pi = pi > y/2Vo and p2 the momentum of this trajectory before 
and inside the barrier respectively, energy conservation gives pl/2 = P2/2 + Vq. This is the 
first equation connecting pi to p2, but we need a second one which is imposed by the 
propagation time T = ti + t2 where: 

ti = is the time to go from q to —a with momentum pi, 

_ xf+a_ ^j^^ time to go from —a to Xj with momentum p2. 
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The combination of these two equations gives 

I + q 



T=-^ + ^I±^ (19) 



which can be rewritten as 

(pj - 2Vo)ipiT + a + qf = (xf + a)Y,. (20) 

This is a quartic polynomial, which we solve numerically. We obtain four solutions: one is 
always negative, which we discard since we fixed the initial position q on the left side of the 
barrier; two are sometimes complex and, when real, have pi < 1; finally, one of the roots 
is always real, larger than 1 and tends to ^^^jA when Vq is negligible (the limit of a free 
particle). We take this last root as the initial momentum pi. 
The action S is also a function of pi given by 



= f t,+ (|-ro)t2 (21) 
1/ \ 1/ \ n? 7^ Vo(xf + a) 

We calculate the derivatives of 5* numerically by computing pi and S for the initial conditions 

dS 

{q,Xf), (g + dq, Xf), {q,Xf + dxf) . . . and approximate — — {z,Xf, T) by [S{z^Xf + dxf, T) — 
S{z,Xf,T)]/dxf, etc. Finally, the propagator inside the barrier is given by Eq. flTT]) plus the 
phase correction C, calculated in the previous subsection. The probability density, which in 
independent of ^, becomes 



(z,x/,T),c|^ = r^^^=^exp 



2 / \ 2 



P-Pl 



(22) 



Figure [3] shows iV'scP as a function of x/ for the same parameters as in subsection IIII A[ 
Although the semiclassical approximation also improves for small h, here we shall fix ^ = 1. 
This is because the behavior of the propagator becomes trivial for small h: if p < 1 the 
wavepacket bounces off the barrier almost completely, and otherwise it simply passes over 
the barrier barely noticing the presence of the potential. 

The first remark is that the wavepacket is continuous at xj = —a: the extra phase 6 
does play its role correctly. As in the case before the barrier, the comparison between exact 
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and semiclassical calculations is always at least qualitatively good, and sometimes even 
quantitatively so. However, there are two main effects that the semiclassical approximation 
cannot take into account. 

1. there is a gap between the exact and semiclassical curves, which decreases progressively 
as Xf increases, and is bridged near the local maximum of the probability density. 
The reason may come from the fact it is not possible to impose the continuity of the 
derivative of ipse with respect to Xf at —a. 

2. there are oscillations on the exact curve (especially for p = 2 and T = 50) close to 
the right side of the barrier, that are not present in the semiclassical approximation. 
This is a purely quantum effect, because classical mechanics can not account for a 
reflected trajectory which would interfere with the direct one in this case, {ipsd"^ is 
in fact the mean-value of the oscillations, and that is why there is a discontinuity of 
the wavepacket at Xf = a, since the exact curve is beginning at the bottom of an 
oscillation. 

If we want to stay strictly in the semiclassical limit, there is nothing we can do about the 
lack of interferences in the barrier region: this is the limit of our approximation. But if we 
want to use the semiclassical point of view in order to provide a more intuitive picture of 
the quantum world, we can add a 'ghost' trajectory that reflects at x/ = a and see if it can 
account for the interferences. Similar ideas have been applied to the frequency spectrum of 
microwave cavities with sharp dielectric interfaces and, more recently, to the spectrum of 
step potentials confined by hard walls lisj] . The argument will be the same as in subsection 
IIII At except of course that the reflected trajectory will now bounce on the right side of the 
barrier. The equation for pi = pi is again a quartic polynomial given by 



We know that pi direct is the same as Pi reflected at Xf = a and we choose the only solution of 
( l23l) which satifies this condition. The expression of the new action is: 



The expressions of ipd and ipr are the same as eq. (fTT!) but with pi, S, nigg and rrigp indexed 
by d or r. After some calculations, the new expression of the probability density inside the 




(23) 




(24) 
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barrier becomes 



{z,Xf,T), 
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exp 
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(25) 
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X exp 
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2 mL J + 



P~Pid 
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P-Pli 



"qqd ' ' qpd 

where 9' is the new extra phase (that absorbs the previously computed ^) and 



rriqpd 



Sf — Sd 1 

<Pr - <^d = T h - arctan 

h 2 \mq„d 
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1 mqqdmqp d fp-pid 

2 ■ + ml 
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2 J, + ^ 



c 



(26) 



qqd ' ■ ■ "gp f 

The results of such an expression, however, are not good: the oscillations become too big, 
which means that the reflected trajectory needs to be attenuated by a reflection coefficient 
p. To calculate p we use the following reasoning: for each point Xf inside the barrier there 
corresponds a reflected trajectory from qio xj with a certain value of pi > 1 computed with 
Eq. (1251) . We take for p the same attenuation coefficient a plane wave with momentum pi 
would have. Let (^F e^'^^f + Ge"***^^) and C e^^^^ be such a plane wave inside and after the 



barrier respectively, where k = ^/2{E — Vo)/h = \/]y{ — p^ jfi and k = V2E/h = pi/h. The 
continuity of this function and its derivative at xj = a give us the relative weight of the 
reflected trajectory with respect to the direct one: 

1-K/k 1 - ^1 - pypj 



piPi,Vo) 



(27) 



1 + K/k 1 + y/l-pypj 

The expression for the total propagator becomes [ipscd + pi'scr^^^')- We use the same 



argument as in subsection IIII Al to compute the extra phase 6', adding another correction 
to the wavefunction on the right side of the barrier. Because there is always a single 

trajectory on the right side, ^' does not affect the probability density there. We find that 

cos^' = Al/2 - 1 where A+ = W{a)/D{a). 

The new results are displayed in figure HI The gap is still present, but the agreement 

between exact and semiclassical on the right side is nearly perfect. The interferences are 
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indeed coming from a real 'ghost' trajectory that bounces off at the end of the barrier 
hke a quantum plane wave. Since the left side of the figure has not changed much, the 
reflected trajectory has no effect on this part of the wavepacket and we don't need to 
consider additional reflections. Furthermore, we don't have to take pipscr into account when 
we calculate 6 in subsection IIIl A[ We finish this subsection with two comments: first, the 
approximation with the ghost trajectory is accurate even for larger values of h. Second, the 
wavepacket becomes continuous at x/ = a. That is very interesting because continuity comes 
only when we include the reflected trajectory, whereas the part of the wavepacket which 
goes through the barrier is calculated independently with a single direct trajectory (see next 
subsection) . This means that the semiclassical propagator after the barrier somehow knows 
there is a reflected part. 

In the next subsection, we will briefly present the computation of the wavefunction at 
the right side of the barrier. 

C. After the barrier: a < x 

Following the same arguments as in subsection IIII B\ we use the energy conservation 
p1/2 = P2/2 + Vq = pI/2 (the index 3 refers to the right of the barrier) and the different 
times t-i = — f = ^ and to = to calculate the initial momentum of the direct 
trajectory. We obtain 



In this region, no reflection is possible and the probability density iV'scP is simply given by 
Eq. (1221) . The results are presented in figure [5l For any values of p, T or h, there is still 
a very small difference between the exact and semiclassical curves for the ascending part of 
the wavepacket, whereas the agreement is perfect when the function is decreasing. 

The conclusion of this section is that the semiclassical approximation with real trajectories 
gives very good results and is indeed able to describe some important quantum effects. 




(28) 



whereas the action becomes 




(29) 
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Interference on the left side of the barrier appears naturally when the wavepacket hits the 
barrier and the comparison with the exact solution gets better as h gets smaller. However, 
these interferences cannot be obtained in the barrier region, since there are no reflected 
trajectories in the classical dynamics. We showed that these interferences can be recovered 
if a 'ghost' trajectory that reflects at x = a is added and assumed to contribute with the 
same coefficient of a plane wave of initial momentum pi. With this addition the semi classical 
approximation becomes again very accurate inside the barrier. In the next section we shall 
briefly discuss the possibility of using our results to calculate the tunneling time as defined 




IV. SEMICLASSICAL TUNNELING TIMES 



The question of how much time a particle spends in the classically forbidden region 
du ring the tunneling process has been attracting the attention of physicists for a long time 



i la, llTl, uM, uM, 12Q|, l2l|, l22|. The very concept of a 'tunneling time' is, however, debatable 



181]. Nevertheless, in a semiclassical formulation where real trajectories play crucial roles in 



the tunneling process, the temptation to estimate such a time is irresistible. 

Since we are considering a wavepacket, and not a classical state localized by a point in 
the phase space, we can only define a mean value of the tunneling time. Let us fix the initial 
conditions q,p (such that p < 1) and Xf > a. The probability of finding the initial Gaussian 
state at X/ after a time T is given by | < Xf\K{T)\z > p. Therefore, the particle can 
reach xj from {q,p) in several different time intervals T. For each value of the time T there 



corresponds a single real trajectory whose initial momentum Pi{T) > p = a/2Vo is given by 
Eq. ( l28l) . This trajectory spends a time t(T) = ^^^^ = ^ in the region —a < x < a. 

Notice that the average energy of the wavepacket is below the barrier but the contributing 
trajectory always has energy above the barrier. Therefore, for fixed q,p,Xf, the probability 
that the wavepacket crosses the barrier in a time t(T) is proportional to | < Xf\K{T)\z > p. 
Following ref. j^, we can define the mean value of the tunneling time as 

r+oo 

{t) = ^-^ T{T)\<x\K{T)\z>\'^dT (30) 
Jo 

where 

r+oo 

m= I < x\K{T)\z > l^dT (31) 

Jo 
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is the normalization factor. It is not equal to 1 because only the part of the wavepacket which 
goes through the barrier is considered. This is important in our case, since the semiclassical 
approximation is better for x/ > a. 

We calculated these integrals numerically, performing a discrete sum over T„ = n5T , with 
n = 1, 2, . . . , and 5T = T^ax/N. If an observer is placed at a fixed position Xf > a, as 
the time T slips by, he/she sees the wavepacket arriving from the barrier, becoming bigger 
and bigger, reaching a maximum and then decreasing and disappearing. We ended the sum 

at Tmax such that | < x\K{T)\z > P < 10"^ V T > Tmax- 

An important remark is that (r) is independent of the observer's position Xf (except for 
small fluctuations due to the numerical computation), since Eq.f l3UI) measures only the time 
inside the barrier. The three different times we are going to use for comparison are: 

• {Tbarr) IS the tunneliug time computed according to Eq. ( l30l) 

• {'Tfree) IS obtained from the same way as {rtarr), but in a system without barrier; (rfree) 
is simply the time for a free wavepacket to go from —a to a. 

• Tc/ass = / ^" is the time required by a classical particle to cross the barrier. 

Fig. [6] shows the dependence of these functions with respect to p. The curves become 
very similar as p increases, because the barrier becomes more and more negligible. The 
wavepacket spreads but stays centered around p, which explains why it behaves like a particle 
of momentum p. When the influence of the barrier is more important, the wavepacket gets 
trapped by the barrier and slows down ((r^arr) is above (Tfree)), but for p < 1.8, (r^ree) and 
Tdass start to increase very fast ( r^ass actually diverges at p = 1), whereas {rbarr) stays 
finite until p is very close to 0: thanks to the tunnel effect the wavepacket is accelerated by 
the barrier, which acts like a filter for the wavepacket and cuts off the contributions of its 
slowest components (see Fig. [71(a)). On the other hand when p increases, the fraction of 
the trajectories with p < V2Vo = 1 becomes negligible and the barrier simply restrains the 
propagation of the wavepacket (fig. [71(b)). 

V. CONCLUSION 

In this paper we used the semiclassical approximation Eq. (ITTl) . derived in j^, to study 
the propagation of a wavepacket through a finite square potential barrier. One of the main 
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purposes of this work was to test the vahdity and accuracy of the approximation, which 
involves only real trajectories, in the description of tunneling. Surprisingly, we have shown 
that the approximation is very good to describe the wavepacket after the barrier, even 
when the average energy of the wavepacket is below the barrier height. The region before 
the barrier is also well described by the approximation, although discontinuities are always 
observed because of the sudden disappearance of the reflected trajectory. The continuity of 
the wavefunction between this region and the region inside the barrier also depends on the 
calculation of an extra phase 6. Finally, inside the barrier the semiclassical formula is not 
able to describe interferences. These, however, can be recovered when a ghost trajectory, 
that reflects on the right side of the barrier, is included and attenuated with the proper 
coefficient. In all regions the approximation becomes more accurate as h becomes smaller. 

The semiclassical approximation ffTTl) is particularly relevant because the propagated 
wavepacket is not constrained to remain Gaussian at all times, as in the case of Heller's 
Thawed Gaussian Approximation 121], and also because it uses only a small number of real 
trajectories. These are much easier to calculate than complex ones, especially in multi- 
dimensional problems. The demonstration of its ability to describe tunneling and interfer- 
ences is important to establish its generality and also to provide a more intuitive under- 
standing the processes themselves. In particular, using the underlying classical picture, we 
have computed a tunneling time which shows that the wavepacket can be accelerated or 
restrained by the barrier depending on the value of the initial central momentum p. 

Some interesting perspectives of this semiclassical theory are the study of propagations 
through smooth potential barriers (which are more realist and more adapted to semiclassical 
calculations), the study of time dependent barriers and the extension of the method to 
higher dimensions and to chaotic systems. 
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FIG. 1: Direct and reflected trajectories from q to xi < —a. For —a<X2<aoTXs>a only the 
direct trajectory exists. 
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FIG. 2: (Color online) Exact (blue thin lines) and semiclassical (red thick lines) wavepacket at 
time T = 50, except for panel (h) where T = 10. We fixed h = 1 for (a), (b) and (c), whereas 
h=l/4: for (d), (e), (f), (h) and h = 1/10 for (g). 
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FIG. 3: (Color online) Exact (blue thin lines) and semiclassical (red thick lines) wavepacket inside 
the barrier for various values of p and T and h= \. 
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FIG. 4: (Color online) Exact (blue thin lines) and semiclassical with ghost reflected trajectory (red 
thick lines) wavepacket inside the barrier . The panels on the right are magnifications of the left 
ones, showing the perfect match between the approximation and the exact solution. 
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FIG. 5: (Color online) Exact (blue thin lines) and semiclassical (red thick lines) wavepacket after 
the barrier for p = 2 and T = 50 and 100. 
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FIG. 6: (Color online) Tunneling time as a function of p. The red thick line and the green curve 
show the semiclassical result according to Eq. ()30p for the square barrier and the free particle 
respectively. The thin black line is the classical time for the square barrier potential. 
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FIG. 7: (Color online) Exact (blue thin lines) and semiclassical (red thick lines) wavepacket after 
going through the barrier. The green curve (above the other two in both figures) shows the 
corresponding free particle wavepacket. For p = 0.5 the barrier acts like a filter and only the fast 
components of the initial wavepacket go through. For p = 2, on the other hand, the wavepacket 
interacting with the barrier is slightly behind the free particle wavepacket, showing that the barrier 



slows the trajectories because the momentum is reduced to p2 = y/p'f — 2Vo < pi between —a and 
a. 
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